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RADIO POSITIONING SYSTEM 



The present invention relates to a radio positioning system. 

EP-A-0 303 371. the contents of which are hereby incorporated by reference, describes 
a radio navigation and tracking system which makes use of independent radio 
transmitters set up for other purposes. The signals from each transmitter, taken 
individually, are received by two receiving stations, one at a fixed and known location, 
and the other mounted on the mobile object whose position is to be determined. A 
representation of the signals received at one receiving station is sent via a link to a 
processor at the other receiving station, where the received signals are compared to 
find their phase differences or time delays. Three such measurements, made on three 
widely spaced independent transmitters, are sufficient to determine the position of the 
mobile receiver in two dimensions, i.e. its position on the ground. The phase or time 
offset between the master oscillators in the two receivers is also determined. 

"CURSOR", as the system described in EP-A-0 303 371 is known, is a radio 
positioning system which can use the signals radiated by existing non-synchronised 
radio transmitters to locate the position of a portable receiver. Unlike some other 
systems which use the temporal coherence properties of networks of purpose-built 
synchronised transmitters, CURSOR makes use of the spatial coherence of the signals 
transmitted by single transmitters. In a further development (see EP-A-0 880 712 & 
WO-A-99/21028), the technology has been applied to find the position of a mobile 
phone handset in a GSM or other digital telephone system, and these are examples of 
an 'Enhanced Observed Time Difference' (EOTD) method using the down-link 
signals radiated by the network of Base Transceiver Stations (BTS) of the telephone 
system. 

In the digital mobile telephone application described in EP-A-0 880 712, the signals 
from each BTS within range of the handset are received both by the handset itself and 
by a fixed nearby receiver, the Location Measurement Unit (LMU), whose position is 
accurately known. Representations of the received signals are passed to a Mobile 
Location Centre (MLC) where they are compared in order to find the time difference 
between them. In Figure 1 we show the geometry of a standard two-dimensional 
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system. The origin of Cartesian co-ordinates x and y is centred on the LMU positioned 
at O. The orientation of the axes is immaterial, but may conveniently be set so that the 
y axis lies along the north-south local map grid. The handset, R, is at vector position r 
with respect to the LMU position O. A BTS, A, is shown at vector position a. 

Consider first the signals from BTS A. The time difference, A/ a , measured between the 
signals received at O and R is given by 

A/ a =(|r-a|-|a|)/u + e, 

where v is the speed of the radio waves, and e is the clock time offset between the 
clocks in the receivers at O and R. The value of € represents the synchronisation error 
between the measurements made by the two receivers. Similarly, we may write for two 
other BTSs (B and C) at vector positions b and c (not shown): 

A/ b =(|r-b|-|b|)/u + c, 

A/ c =(|r - c| - \c\)/u + e . (1) 

The values of A/ a , At b , At c , are measured by the methods disclosed in EP-A-0 880 712 
and the values of a, b, c, and u are known. Hence the equations (1) can be solved to 
find the position of the handset, r. 

In WO-A-99/21028 we describe how these same time offsets can be measured using 
locally-created templates in a GSM telephone system as follows. Suppose that the 
handset has recorded a short burst of the GSM signals from BTS A. Contained within 
that recording is the framing structure, synchronisation bursts and other 'given' data 
(or predetermined values) which are a constant feature of those transmissions. The 
processor within the handset can create a matching template, based on the known 
structure of the network signals. Received signals can then be matched by the locally- 
generated template. When the template finds a match, the correlation peak at the 
position of best match corresponds to the time offset between the received signals and 
the local clock inside the handset. This time offset, Ar aJ is given by 



A/ai =(|r - a|)/u + a a + e ] . 

where a a is the time offset of the BTS transmissions and £j is the time offset of the 
handset's internal clock, both relative to an imaginary universal 'absolute' clock. The 
signals from BTSs B and C may also be measured in the same way, giving 

A/ b] =(|r - b|)/u + a b + e u 

and 

A/ C ] =(|r - c|)/u + a c + e x . (2) 
The same measurements can also be made by the LMU, giving 
A/a2 =(N)/o + a a + e 2 , 
A/ b 2 =(|b|)/u + a b + e 2 , 

and 

A/ C 2 =(|cl)/u + a c + € 2 , (3) 

where e 2 is the time offset of the LMU's internal clock relative to the imaginary 
universal absolute clock. Subtracting equations 3 from 2 gives us 

A/ a = A/ a i - A/a2= (|r - a| - |a|)/u + e , 

A/ b = A/ b i - A/ b2 = (|r - b| - |b|)/u + £ , 

and 



A/ c = A/ c , - At c2 = (|r - c| - |c|)/u + £ , 



(4) 



where € = €\ - s 2 . It will be noted that equations 4 are just like equations L and can be 
solved in the same way to find the position of the handset, r. 

It will be apparent that the CURSOR method, in common with all other EOTD 
methods, requires a network of LMUs to be set up within the coverage area of the 
telephone system. These units act as reference points at which the unsynchronised 
signals radiated by the BTSs are measured for comparison with the same signals 
received by a handset. Each position measurement requires a match to be made 
between the signals received by the handset from a number of nearby BTSs, and 
signals received by an LMU from the same set of BTSs. In practice, it is often difficult 
to find a match using just one LMU, especially if the LMU network is sparse. It is 
therefore necessary to combine the measurements from two or more LMUs. However, 
each new LMU brought into the calculation adds a further unknown clock time offset 
(£ 2 , £3 etc.), each of which therefore requires an additional BTS to be measured to 
provide the extra equations needed to solve for all the unknown quantities. A 
minimum set of measurements for positioning in two dimensions requires at least three 
geometrically-dispersed BTSs in the set if one LMU is involved. Adding a second 
LMU increases the minimum number of BTSs required to four, and so on. 

One solution to this problem is presented in WO-A-99/21028 where it is shown how 
the LMU network can be synchronised. Suppose that an adjacent pair of LMUs, Lj 
and L 2 , can see a common BTS. The positions of the LMUs and the BTS are all 
known, so a single measurement of the BTS signals by each LMU is sufficient to 
determine the clock time offset between the LMUs. For example, suppose that the 
distance from L\ to the BTS is s\ 9 and the distance from L 2 to the BTS is s 2 . L] 
measures time offset At] and L 2 measures A/ 2 , given by 

Ari = S]/v + a + €2 , 



A/ 2 = s 2 /u + a + £ 3 , 



where a is the time offset of the BTS transmissions, and e 2 and £ 3 are the time offsets 
of the LMU internal clocks. Subtracting the second equation from the first yields 

e 2 - €3 = At) - At 2 + s\/v - S2/0 5 

which is the relative time offset of Li with respect to L 2 . This process may be repeated 
for a second pair of LMUs, say L2 and L3, and another BTS whose signals can be 
received by both of these second pair of LMUs. In this way a synchronisation map 
may be calculated, which provides the clock offsets of all the LMU internal clocks 
relative to one of them adopted as the master 'LMU network clock time'. Having 
established the LMU synchronisation map in this fashion, a CURSOR position 
measurement can include any number of LMUs without the penalty of adding an extra 
unknown quantity for every LMU. 

The present invention teaches how the same advantages of a synchronised LMU 
network may be obtained by setting up one or more 'virtual LMUs 5 in the network 
which act as interface nodes for the real LMUs. 

According to a first aspect of the invention, there is provided a method of generating a 
list of the timing or phase offsets of a plurality of transmission sources relative to a 
common timing reference, the method comprising 

(a) acquiring data from one or more receivers, the positions of which are 
known or determined, the data from a receiver comprising timing or 
phase offsets of signals received from the transmission sources relative 
to a timing reference source or to each other; and 

(b) combining the acquired data and producing a list of the offsets, relative 
to the common timing reference, of a number of the transmission 
sources. 

In practice, particularly when the transmission sources are transmitters in a digital 
mobile telephone network, the offsets from the list are used in place of offsets 
obtained directly from the receiver or receivers. 



In a method using techniques similar to or as described in EP-A-0 880 712 ; in place of 
timing or phase offsets, data representative of the received signals from which the 
timing or phase offsets of signals received from the transmission sources relative to 
the timing reference source may be used. 

5 

Therefore an alternative first aspect of the invention comprises a method of generating 
a list of the timing or phase offsets of a plurality of transmission sources relative to a 
common timing reference, the method comprising 

(c) acquiring data from one or more receivers, the positions of which are 
10 known or determined, the data from a receiver being representative of 

the received signals from which the timing or phase offsets of signals 
received from the transmission sources relative to a timing reference 
source or to each other may be determined; and 

(d) combining the acquired data and producing a list of the offsets, relative 
15 to the common timing reference, of a number of the transmission 

sources. 




The invention also includes apparatus (a 'virtual LMU') for carrying out either or both 
of these methods. The apparatus may include a computer located at any one of a 
20 number of places connected to the network and programmed to carry out the required 
process. 



A virtual LMU is a computer process which may run on any computer platform 
connected to the network. It is assumed that the network of BTSs is unsynchronised in 

25 that the relative phases (the relative transmission delays) of the BTS signals bear no 
constant or known relationship to each other, but that nevertheless the BTS oscillators 
are quite stable, so that their relative phases change only slowly with time. In these 
circumstances, it is possible to predict the current relative phases from sufficiently- 
recent historical data. The real LMUs in the network make measurements of all the 

30 BTSs they can detect in a cyclic fashion, repeating the cycle every few seconds. They 
maintain these measurements in a stack, replacing the oldest measurements with the 
most recent. A linear or low-order polynomial fit to the measurements therefore 
provides a predictor for extrapolation into the near future, or for interpolation in the 
recent past. Let us suppose that the phases are sufficiently stable that reliable 
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predictions can be made over a period of say, ten minutes. Then every ten minutes, 
the virtual LMU (VLMU) contacts each real LMU and receives its predictors for the 
BTSs in its measurement set. It is likely that many of the BTSs will have been 
measured by more than one LMU, so the VLMU analyses the complete data set from 
all the real LMUs to determine the best values of the LMU clock time offsets and the 
BTS relative phases using well-known methods. 

By way of example, one such method is now described. Consider a network of N 
LMUs and M BTSs. Let the position of the n th LMU be represented by the vector L n 
and let the position of the m th BTS be represented by the vector B m , both measured 
with respect to the same origin. Signals radiated by BTS m will be received by LMU n 
after a time lag, A/ nm5 where 

A/ nm = ILn-BJ/o + € n + a m ± CTnm , 

£„ is the clock offset of LMU n, and is the clock offset of BTS m, both with respect 
to an imaginary universal 'absolute 1 clock, and cr nm is an estimate of the error in the 
measurement of Ar nm . Over the entire network of N LMUs, all of the M BTSs are 
visible. Each individual LMU, however, will only see a few of them, but as long as 
there is significant overlap of visibility, it is possible to take the set of all At values and 
solve for values of ^ and a m . Hence the VLMU can calculate timings for any BTS as 
if the network of LMUs were synchronised. 

To illustrate this, a simplified problem is shown and solved below using N = 2 and M 
= 4. We choose = 0 for simplicity. We are allowed to do this as the 'absolute' clock 
time is completely arbitrary and may, for example, be measured by the internal clock 
of LMU number one. The equations can be written in matrix form 
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or equivalently as A.x = b + V, where V is an unknown vector of the actual errors on 
each measurement. Let the matrix W be defined by 
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The standard technique known as 'least squares' postulates that the estimate for x 
which minimises V is given by 



x = (A T WA)" 1 A T Wb . 



This particular example can be solved explicitly. For simplicity, let us assume that all 
the values of a nm are the same, and equal to o. This gives the result 



a, =(3Ai +A2 + -D21 
a 2 = (3A 2 + Ai +£>22 

cx 4 = (2D 2 4 + Dn+D n 



D22 )/4 ± 0.87a, 

D 2i )/4 ±0.87ct, 
± l.OOo-, 

D 2i -D22)/2 ±\Ala, 
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£2= (An + D22 - D } ] - D }2 )/2 



± 1.00a, 



L - B 

where £> = At - - 

O 



Note that even in this simple case, when a BTS is seen by both LMUs 5 the errors in the 
calculated offsets are less than those in each of the measurements themselves. This is 
an important advantage of the virtual LMU method. 

The LMUs may also contain other means of synchronisation. For example, each real 
LMU could be connected to a GPS or other timing reference receiver which serves to 
provide the common timing reference. In this case, the LMU network may be 
considered as synchronised already to this common timing reference (say GPS 
standard time), and then the VLMU need not solve for the individual values of e as 
these are already known. 

The VLMU mode of operation described above may be called the 'pull mode' as it 
requires the VLMU to instigate data transfer to itself from every real LMU. It is also 
possible to have each real LMU continuously check the difference between its own 
prediction of every BTS phase using the last predictor sent to the VLMU and the 
actual measured phase. When this difference exceeds a given value, it can send its new 
predictor set to the VLMU. This mode of operation may be called the 'push mode 9 . 
The particular mode appropriate for a real system depends, amongst other things, on 
the stability of the BTS network. 

A second aspect of the present invention shows how the CURSOR method (or other 
EOTD method) can be applied without the need of a network of real LMUs. 



According to the second aspect of the invention there is provided a method of 
determining the position or change in position of a receiver or receivers in a network 
of transmission sources, the method comprising 



(a) at a first time, measuring the relative time or phase offsets with respect 
to each other or with respect to a timing reference source of the signals 
received by the receiver from a plurality of the transmission sources; 

(b) at a selected second later time, measuring the relative time or phase 
offsets with respect to each other or with respect to a timing reference 
source of the signals received by the same or another of the receivers 
from a plurality of the transmission sources; 

(c) calculating the relationship between first and second curves 
representing the loci of possible positions of the receiver or receivers at 
the first and second times respectively; and 

(d) determining the position or change in position of the receiver or 
receivers at the first or second times from the positions of the 
intersections of the curves. 

In the method defined above, the receivers may be the same receiver. 

The second aspect of the invention also includes a method of determining the position 
or change in position of one of a plurality of receivers, the positions of all of which are 
unknown, in a network of transmission sources, the method comprising 

(e) at a first time, measuring the relative time or phase offsets with respect 
to each other or with respect to a timing reference source of the signals 
received by the receiver from a plurality of the transmission sources; 

(f) at a selected second time, measuring the relative time or phase offsets 
with respect to each other or with respect to a timing reference source 
of the signals received by the same or another of the receivers from a 
plurality of the transmission sources; 

(g) calculating the relationship between first and second curves 
representing the loci of possible positions of the receiver or receivers at 
the first and second times respectively; and 

(h) determining the position or change in position of the receiver or 
receivers at the first or second times from the positions of the 
intersections of the curves. 

In the method defined above, the first and second times may be the same time. 
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The invention also includes apparatus for carrying out either or both of these methods, 
which apparatus may include receivers suitably programmed. 

To illustrate one example, in which the transmitters and receivers are part of a digital 
mobile telephone network, let us suppose that the handset is at vector position r(/]) at 
time /]. Equations (2) are then 

A' a] (/i) = (|r(/,) - a|)/u + a a (r,)+ £,(/,) , 

A/bi('i) = (|r(/0 - b|)/u + a„(/,)+ e^u), 

Af c ,(r,) = (|r(/i) - c|)/t» + Oc(/,)+ £,(/,), (5) 

where £ a 0i) is the time offset of the BTS transmissions from A, £b('i) the offset from 
B, £ c (/i) the offset from C, and e^ti) is the time offset of the handset's internal clock, 
all measured at time /] relative to an imaginary universal 'absolute 1 clock. Now let the 
handset be at a different vector position r(t 2 ) at a later time / 2 . Equations (2) are then 

A/ a i(/ 2 ) = (|r(/ 2 ) - a|)/u + Oa(/ 2 )+ £i(/ 2 ) , 

A/bifo) = (|r(/ 2 ) - b|)/u + a b (r 2 )+ e,(/ 2 ), 

Ar c ,(/ 2 ) = (|r(/ 2 ) - c|)/u + ac(r 2 )+ e,(/ 2 ). (6) 
Subtracting equations (6) from equations (5) gives 

A/ai(/i) - Afaifo) = (|r(/,) - a| - |r(/ 2 ) - a|)/u + (a a (/,) - Oa(r 2 ))+ (£,(/,) - £l (/ 2 )) , 
A/bi('i) - A/ bl (/ 2 ) = (|r(/,) - b| - |r(Y 2 ) - b|)/u + (a„(/,) - a b (/ 2 ))+ - s,(f 2 )), 

A/ c i(/i) - A/ c ,(/ 2 ) = (|r(/,) - c| - |r(/ 2 ) - c|)/u + (««.(/,) - Oc(/ 2 ))+ (£,(/,) - e,(/ 2 )). (7) 
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As has been noted above, the BTS transmissions in a GSM or other digital network are 
not usually synchronised. However, they are often derived from a common reference 
timing source. Where this is the case, the signals therefore keep a constant phase offset 
with respect to each other, so e^t } ) = s a (t 2 ), £bOi) = s b {t 2 \ and £ c (t } ) = £ c (/ 2 ). (Even 
where this is not the case, the BTS signals are often sufficiently stable to allow an 
estimate to be made of these differences in a VLMU or other network element.) 
Writing A/ a = A/ a] (/,) - Ar al (/ 2 ), At b = A/ b i('i) - Af bl (f 2 ), A/ c = A/ Cl <7,) - A/ c i(/ 2 ), and € = 
- £](t 2 ), we get 

Af. = (|r(r,) - a| - \rit 2 ) - a|)/u + e, 

At h = (|r(/0 - b| - |r(/ 2 ) - b|)/u + e, 

At c = (|r(/,) - c| - |r(/ 2 ) - c|)/u +e. (8) 

Equations (8) can be understood graphically by reference to Figure 2. The diagram 
represents the plane of the Earth's surface (assumed to be flat) near to BTS units A 
and B. In the first of the equations, the term |r(/j) - a| represents the distance of the 
handset from A at time tu and the term |r(/ 2 ) - a| represents its distance at time t 2 . The 
first of equations (8) can be rewritten as 

|r(/j) - a| - |r(/ 2 ) - a| = uA/ a - ve, 

which therefore represents the loci of two concentric circles centred on A whose radii 
differ by uA/ a - vs. These are marked 1 and 2 respectively in Figure 2. The above 
equation does not define the radius of either circle, but only the difference between 
them. By itself, this equation therefore does nothing to locate either of the points r(t\) 
or r(r 2 ). Also marked in Figure 2 are a second pair of concentric circles, 3 and 4, 
representing the second of the equations (8) and centred on B. Again, their radii are 
unconfined by the equation, but the difference between their radii must be vAt b - vs. 
Point r(/i) must lie at one of the intersections of circles 1 and 3, and point rfo) must 
lie at one of the intersections of circles 2 and 4. Let us suppose that the value of e is 
zero, i.e. that the handset's internal clock has kept perfect time between the 
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measurements. Then if we know the position of point r(/j), say P in Figure 2, we can 
deduce that point r(/ 2 ) must be at Q, because the four circles are now fixed in space, 
and hence we have measured the position of the handset at the later time. In practice, 
we can't assume that e is zero, so we must use all three of equations (8) to find Q 
given P. 

It will be noted there is also a second point of intersection, Q' ? of the circles 2 and 4, 
and hence an ambiguity in the determination of r(f 2 ). If this ambiguity cannot be 
resolved by other means (for example, by knowing that the handset is being used by a 
pedestrian who could not have moved from P to Q* in time t 2 - t\) then four 
measurements involving four BTSs must be used. 

One of the advantages of the second aspect of the present invention is that if a single 
handset makes measurements of the signals from a minimum of three geographically 
separate BTSs at two different times, then it is possible to determine the change in the 
position of the handset between these two times (without reference to a known starting 
point), provided that at least three of the BTSs are common to the two sets of 
measurements. If the position of the handset has been determined at some point 
previously, use of this technique may enable applications requiring autonomous 
navigation within the handset, i.e. navigation without further reference to the 
CURSOR or any other EOTD system. 

The above discussion shows how the position of a handset at a later time can be 
obtained from three CURSOR-like measurements on three geographically-dispersed 
BTSs made both at the later time and at an earlier time provided (a) that the position of 
the handset is known at the earlier time, and (b) that there is no relative drift between 
the signals radiated by the BTSs in the interval (or that such drift is known). The 
accuracy of the position determination depends both on the precision with which the 
first position is known and on the distance moved between the measurements. 
Equations (8) actually contain five unknown quantities: the x and y co-ordinates of 
each of the two points P and Q, and the unknown handset clock drift £. CURSOR-like 
measurements on each of five geographically-dispersed BTSs at both a first and a 
second position are therefore sufficient to determine both P and Q uniquely if the 



handset has moved appreciably between the two sets of measurements. Here, then, is 
an EOTD method of finding the position of a moving handset in an unsynchronised 
digital mobile phone network without the need of any LMUs at all. It may be 
especially useful for tracking the handset in a network which does not already have a 
location system. 

A further advantage of the second aspect of the present invention is to use the known 
position of one handset to find the unknown position of a second handset. Suppose 
that each handset makes a measurement at about the same time using three 
geographically-dispersed BTSs. Equations 5 for the first handset, Hi, are then 

A/aiCO = (InCi) - a|)/u + a a (r,)+ <?,(/,) , 

AfbiCi) = (|ri(fi) - b|)/u + Ob(/i)+ ei(ti), 

A/ciCi) = - c|)/u + Oc(/,)+ *,(*,), (9) 

where rj(/i) is the (known) position of H] and the measurements are made at time t\. 
The same equations for the second handset, H 2 , are 

A/ a i(/ 2 ) = (|r 2 (/ 2 ) - a|)/u + a a (/ 2 )+ e 2 (t 2 ) , 

A/ b ite) = Qr 2 (t 2 ) - b\)/o + oe b (/ 2 )+ e 2 (t 2 \ 

A/ C i(r 2 ) = (|r 2 (r 2 ) - c|)/u + Oc(/ 2 )+ e 2 (t 2 \ (10) 

where r 2 (t 2 ) is the (unknown) position of H 2 and the measurements are made at time / 2 . 
Subtracting equations (10) from equations (9), noting that a a (/i) = ajj 2 \ a b (/i) = 
a>b(t 2 \ and Oc(fi) = ajj 2 \ and writing A/ a «= A/ a i(/0 - Af a i(/ 2 ), A/ b = A/ bJ (^i) - A/ bI (f 2 ), 
At c = A/ c i(/i) - A/ c i(/ 2 ), and e = s } (t\) - e 2 (t 2 \ we get 

A/. = (|r,(/,) - a| - \r 2 (t 2 ) - a|)/u + 
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A/ b = (ir,(/,)-b|-|r 2 (/ 2 )-b|)/u + c > 



A'c = (|r,(/j) - c| - |r 2 (/ 2 ) - c|)/u +e, (11) 

which are identical to equations (8) in form. Hence, given nfo), the position of Hj at 
time /], we can calculate r 2 (/ 2 ) which is the position of H 2 at time t 2 . (Again, 
measurements on four BTSs are needed to resolve the ambiguity in r 2 (/ 2 ) if this cannot 
be resolved by other means.) In this case, t\ may be equal to t 2 . 

This same idea can be extended to many handsets. In equations (1 1), the vector rjfo) 
refers to a first handset whose position is known, and the vector r 2 (/2) refers to any 
second handset, perhaps representing any one of a large number. It is therefore 
possible to use, as a temporary measure, a first handset at a known location as the 
LMU, offering the possibility of establishing an EOTD service very quickly in a new 
area. 

We can also use this method even in an area already covered by an LMU. Having used 
a normal EOTD method to find the position of a handset, it (the handset) can then be 
tracked without further reference to the LMU. The calculation of the handset's 
position could be made at the MLC, at another location, or in the handset itself, 
depending on the application. 

Another interesting application is to a system of many handsets in an area in which 
none of their positions is known. Provided that any pair sufficiently far apart can 
measure the same five BTSs, their positions can be determined using five equations 
like the three in (8). Otherwise, provided that separated pairs can measure at least three 
or four common BTSs, calculations based on many such measurements by pairs over a 
short period may be sufficient to determine the positions of all of them. Consider, for 
example, a set of three handsets, H l5 H 2 , and H 3 . The relative positions of Hj and H 2 , 
H 2 and H 3 , and H 3 and H], may be computed using three or four BTS measurements 
since the constraints imposed by the equations on the positions of the vertices of the 
triangle HiH 2 H 3 may make the solution unique. 



15 



In the above discussions, we have made the assumption that the BTSs are frequency- 
locked together, i.e. that their signals are all locked to the same primary time standard 
so that their relative transmission delays (RTD) remain constant with respect to each 
other. It is possible to apply the invention to a network in which the BTSs signals are 
5 free-running with respect to each other, i.e. one in which there is no synchronisation of 
any sort between the transmissions. In this case, at least two handsets are required 
which can receive the signals from the same five geometrically-dispersed BTSs at 
substantially the same time as each other. Equations (11) then become 

10 A/. = (|r,(r,) - a| - |r 2 (/,) - a|)/u + €, 

A/ b = (|n(/i) - b| - |r 2 (r,) - b|)/u + e, 

Af c = (|r 3 (/i) - c| - |r 2 (/,) - c|)/u + e, 

15 

At 6 = (|r,(/i) - d| - |r 2 (/i) - d|)/u + e, 

A/ e = (|r,(/,) - e| - |r 2 (/,) - e\)/o +£, (12) 

20 where a, b, c, d, and e are the vector positions of the five BTSs, ri(fi) is the vector 
position of the first handset and r 2 (/i) is the position of the second handset, both sets of 
measurements being made at the same time t\. This moment could be signalled, for 
example, by a particular element of the transmissions such as a particular frame 
number or a special signal. In equations (12), the RTDs cancel out and do not appear 

25 provided that the measurements are made by the two handsets sufficiently close 
together in time that the drifts in the RTDs are small enough to ignore. The five 
measurements in equations (12) are sufficient to locate the two handsets. Note that no 
LMU or equivalent is required in this case. 

30 An extension to this invention is to set up one or more VLMUs in a GSM or other 
transmission network which has no real LMUs at all. Timing measurements by 
handsets like those described above, especially in relation to equations (5) to (12), can 
be processed to provide a map of the BTS relative transmission delays, which may be 
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kept in the VLMUs and used in subsequent standard CURSOR or other EOTD 
position calculations. One might also imagine a system in which a sparse coverage of 
real LMUs is augmented in one or more VLMUs by handset measurements to provide 
the same level of service as from a full network of real LMUs. 
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